Magnetic focusing in normal-superconductor hybrid systems: a semiclassical analysis 
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We study a transverse electron-hole focusing effect in a normal-superconductor system. The 
spectrum of the quasiparticles is calculated both quantum mechanically and in semiclassical ap- 
proximation, showing an excellent agreement. A semiclassical conductance formula is derived which 
takes into account the effect of electron-like as well as hole-like quasiparticles. At low magnetic 
field the semiclassical conductance shows characteristic oscillations due to the Andreev reflection, 
while at higher fields it goes to zero. These findigs are in line with the results of previous quantum 
calculations and with the expectations based on the classical dynamics of the quasiparticles. 
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I. INTRODUCTION 

Investigation of electron-transport properties of 
normal-superconductor (NS) hybrid nanostructures have 
attracted considerable interest both experimentallji 2 ^'^ 
and theoretically^i^ in recent years. In NS hy- 
brid systems a crucial physical phenomenon is the An- 
dreev reflectior>i£, whereby an electron incident on a 
superconductor-normal interface is (partially) retrore- 
flected as a hole into the normal conductor and a Cooper 
pair is created in the superconductor. The first direct 
experimental observation of the peculiar property of the 
Andreev reflection, i.e. that all the velocity components 
are reversed was achieved by Benistant et ali 2 -^ using the 
versatile tool of transverse electron focusing (TEF) 11 . 
The experimental and theoretical investigation of the 
two-dimensional electron gas using the TEF technique 
has been pioneered by van Houten et alH ( see also a 
recent review^ discussing these experiments in terms of 
coherent electron optics). 

Recently, Hoppe et al£ studied theoretically the in- 
terplay of the Andreev reflection and cyclotron motion 
of quasiparticles at the interface of semi-infinite super- 
conductor and normal regions in a strong magnetic field 
parallel with the interface. They found that similarly to 
the normal quantum-Hall systems, edge states are formed 
which propagate along the NS interface but these "An- 
dreev" edge states consist of coherent superposition of 
electron and hole excitations. Therefore they are a new 
type of current carrying edge states which are induced 
by the superconducting pair potential. The authors of 
Ref. i also showed that a semiclassical approximation 
can give a good agreement for the energy dispersion of 
the Andreev edge states with the exact results obtained 
by solving the Bogoliubov-de Gennes equation^ (BdG). 
A clear experimental evidence for the electron and hole 
transport in edge states have been reported by Eroms et 
ali. 

In a disk geometry, it was shown in Refs i 15 ' 16 that such 
edge states can exist both in the presence but also in the 
absence of any magnetic field and that the bound state 



energies calculated semiclassically agree very well with 
the results obtained from the BdG equation. 

Giazotto et al. 8 have extended the study presented in 
Ref. l<5| by considering the effect of the Zeeman splitting 
and of the diamagnetic screening currents in the super- 
conductor on the Andreev edge states. Very recently, Fy- 
tas et ali£ studied the magnetic field dependence of the 
transport through a systems consisting of a normal bil- 
liard and a superconducting island, while Chtchelkatchev 
and Burmistrov have also considered the role of the sur- 
face roughness in NS junctions^. 

In this work we consider the NS hybrid system de- 
picted in Fig. [TJ It is similar to the experimental setup 
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FIG. 1: The hybrid NS nanostructure that we investigate. 
It consists of an infinitely long two dimensional ballistic nor- 
mal conductor of width W coupled to a semi-infinite spin- 
singlet superconductor region. The conductance is measured 
between two normal conducting quantum point contacts: an 
injector (of width Wi) and a collector (of width W c and at 
distance L from the injector). The magnetic field B is ap- 
plied perpendicular to the system (in our calculation B > 
corresponds to a field pointing out of the plane of the system) . 

of Ref. 2 but in our system the normal conducting region 
is a two-dimensional electron gas. It is assumed that the 
quasiparticle transport in the waveguide is ballistic and 
that the waveguide can act as drain which absorb any 
quasiparticles exiting to the left of the injector or to the 
right of the collector. 

First, we calculate the eigenstate of such a system 
where the quantum point contacts are not present, ie 
we consider a normal waveguide in contact with a super- 
conductor. We show how the interplay of the lateral con- 
finement brought about by the finite width of the lead, 
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the applied magnetic field and the proximity effect gives 
rise to a rich physics in this system. We calculate the 
eigenstates of the system by solving the BdG equation 
and then these exact quantum results are compared with 
results obtained from semiclassical calculations. As we 
shall show below the agreement is excellent. We note 
that the semiclassical approximation we used is appli- 
cable in a wider parameter range than the one used in 
Refs. 5.8 for a similar system. 

Having obtained the exact eigenstates and their semi- 
classical approximations, we then turn to the calculation 
of conductance for the case when the two quantum point 
contacts, as depicted in Fig. [TJ are present. The con- 
ductance between the injector and the collector is deter- 
mined by adopting the method used by van Houten et 
alJ£ to describe the TEF in a two dimensional electron 
gas. In our system however, we have to take into ac- 
count the dynamics of both the electron and of the hole 
quasiparticles. Thus, our results on the conductance can 
be considered as a generalisation of the corresponding 
semiclassical calculations of Ref. [l2|. 

Recently, the same NS system has been studied nu- 
merically^ using a Green's function technique^ 8 . The 
influence of the underlying classical dynamics on the con- 
ductance has been, however, discussed only qualitatively. 
Our rigorous semiclassical treatment gives a quantitative 
analysis of the dependence of the transport on the clas- 
sical dynamics of quasiparticles. 

The rest of the paper is organised as follows. In Sec- 
tion [ill we present the exact quantum calculations based 
on the BdG equation. Then in Section HTT1 we discuss the 
results of the semiclassical approximations of the exact 
quantum calculations. 

In Section IIVI we compare the results of the quantum 
and of the semiclassical calculations and give the physical 
interpretation of the results. Section [V] is devoted to the 
semiclassical calculation of the conductance between the 
injector and the collector which is the central result of our 
paper. Finally, in Section [VI] we come to our conclusions. 



II. QUANTUM CALCULATION 

In this section we consider a system consisting of a nor- 
mal conducting waveguide of width W in contact with a 
semi-infinite superconducting region. We derive a secu- 
lar equation solutions of which give the energies of the 
bound states of the system. 

The eigenstates and eigenenergies can be obtained 
from the BdG equation: 
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A* 



(1) 



where ^ is a two-component wave function and Ho = 
(p — eA) 2 /(2m) + V — E-p is the single-electron Hamil- 
tonian (for simplicity, we assume that the effective mass 
m and the Fermi energy Ep is the same in the N and 



S regions^ 9 .). The excitation energy E is measured rel- 
ative to Ep. Scattering at the NS interface is mod- 
eled by an external potential V(x) = Uq6(x). Hard 
wall boundary condition is imposed at the wall of the 
wave guide which is not adjacent to the superconduc- 
tor, i.e. &(x = W,y) = 0. The bound state energies 
are the positive eigenvalues < E < A of the BdG 
equation^. The superconducting pair potential A can 
be approximated by a step function A(r) = A 6(— x) 
without changing the results in any qualitative way 20 . 
Owing to the translational symmetry along the y di- 
rection it is convenient to chose the Landau gauge for 
the vector potential: A(r) = B (0, x, 0) T (here T de- 
notes the transpose of a vector) . Thus, the wave function 
^f(x,y) = (^ c (x,y),^h(x,y)) can be separated and in 
the N region it reads 



* (N) (z,y) 



A e <S>i N \x) 



o ik V 



(2) 



where k is the wave number along the y direction and the 
amplitudes A e will be determined from the boundary 
conditions given below. Substituting ^^(x,y) into the 
BdG equation we find that the function (x) satisfies 
the following one-dimensional Schrodinger equation: 



where 
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(3a) 



(3b) 



Here / = H/\eB\ is the magnetic length, cu c = \eB \ /m is 
the cyclotron frequency, and i>o = 2Ep / (Hlj c ) is the filling 
factor. Equation (|3"a|) is a parabolic cylinder differential 
equation™ and its solution can be expressed in terms of 
the Whittaker functions U(a, £) and V(a,£,): 

= ^(0,0-^4^0, (4) 
V(a,£ w ) 

where £w — V%{W/l — sign(e.B)W). Note that the func- 
tion $i (x) satisfies Dirichlet boundary condition at 
x = W. From the BdG equation it follows that for the 

hole component (x) of the wave function the sym- 
metry relation 



$l N \B,E,x) = & e N H-B,-E,x) 



(5) 



holds. 

The magnetic field is assumed to be screened from 
the superconducting region, hence the vector potential 
is taken to be zero (for the case of finite magnetic pene- 
tration length see e.g. Ref. H). Therefore, in the S region 
the wave function ansatz with eigenenergy E can be writ- 
ten asi£ 



* (S W) = 



c e n°WHx)+cjyH s \x) 



(6a) 



3 



where 



k 2 

q e ,h = k F Jl - p- =f ir}, 



(6b) 
(6c) 
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and fcp 



^/2mE F /U is the Fermi wave number (the upper/lower 
signs in the expressions correspond to the electron/hole 
component). Note that in the S region the wave function 
*&( s \x,y) goes to zero for x — > — oo. 

The four coefficients A e ^, C e .h in Eqs. and (|6"aj) are 
determined by the the boundary conditions at the NS 
interfac o 15 i 16 : 



*(N)| =m(S)\ 

lx=0 \x=0' 



d r 

dx 



*(N) „ ^(S) 



^2 u ° lsc=0' 



(7) 



for all y and /c. The matching conditions shown in Eq. ([7]) 
yield a secular equation for the eigenvalues £ as a func- 
tion of the wave number k. Using the symmetry relations 
between the electronic and hole-like component of the 
BdG eigenspinor given by Eq. {5]), the secular equation 
can be reduced to 15 ' 16 

Im{ 7c D o (E,B)D h (E,B)} = 0, (8a) 

where the 2x2 determinants D e and Dh are given by 



D e (E,B) = 



(S) 



D h (E,B) = D e (—E, 



Z<$> 
-B 



(S) 



(8b) 
(8c) 



Here Z = (2m /ti 2 ) Uq is the normalised barrier strength, 
and the prime denotes the derivative with respect to x. 
All functions are evaluated at the NS interface ie. at 
x = 0. The secular equation derived above is exact in 
the sense that the usual Andreev approximation is not as- 
sumed—. An analogous result was found previousl y 15 ! 16 
for NS disk systems. 



III. SEMICLASSICAL APPROXIMATION 

As a first step to calculate the conductance between 
the injector and collector we should solve the exact 
quantisation condition ([5]) which involves evaluation of 
parabolic cylinder functions. It turns out that for cer- 
tain parameter ranges this makes the actual numerical 
calculations rather difficult. However, as we are going 
to show it in Section IIV1 for the quasiparticle disper- 
sion relations which will be important in the subsequent 
analysis, one can obtain excellent approximations us- 
ing semiclassical methods. The use of the semiclassical 
approximations makes the numerical calculations much 



simpler and gives a better understanding of the under- 
lying physics. The semiclassical calculations are based 
on (i) the WKB approximation^! of the functions $ e (x), 
&h(x) [see Eqs. (I3"a)) and (J5J)] and their derivatives (ii) 
the Andreev-approximation^. The approximated wave 
functions are then substituted into Eq. (I5a|) to obtain 
the semiclassical quantisation conditions. The calcula- 
tions can be carried out in a similar fashion as in Ref. [lrl 
therefore in this section and in the next one we only sum- 
marise the main results. Throughout the rest of the pa- 
per we assume ideal NS interface, ie we set Uq = 0. 

Depending on the energy of the electrons (holes) and 
the applied magnetic field, eight different types of orbits 
can be distinguished. These orbits, denoted by capital 
letters A to H, are shown in Fig. [5] In the geometri- 
cal construction of the classical trajectories we took into 
account that the chiralities of the electron-like and the 
hole-like orbits are preserved when electron-hole conver- 
sion occurs at the NS interface 8 . 

Type A orbits correspond to skipping motion of al- 
ternating electron and hole quasiparticles along the NS 
interface. Neither the electron's nor the hole's trajectory 
hits the wall of the wave guide at x = W. This type of 
orbit was first considered by Hoppe et al£. Type B and 
C orbits are similar to type A but the either the electron 
or the hole can now reach the wall of the wave guide at 
x = W. In case of type D orbits the electrons and holes 
bounce both at the wall of the wave guide and at the NS 
interface, while for type E (G) orbits the quasiparticles 
are moving on cyclotron orbits. Type F (H) is the famil- 
iar edge state of the integer quantum Hall systems. The 
electrons (holes) are moving on skipping orbits along the 
wall of the wave guide at x = W. 

The semiclassical quantisation condition for the orbits 
shown in Fig. [2] can be written in the following simple 
form: 



N(E) = n + /i, 



for n£Z, 



(9a) 



where N(E) can be expressed in terms of the (dimen- 
sionless) classical action integrals S e (r+,Tl), Sh(T+,r!t) 
(see Table H]) for different types of orbits and n is the 
corresponding Maslov index. The actions S e> h are given 
by the following equations: 



S e (rl,Tt) = 2 (e e (r;) - e e (7f )) , 
S h = S e (-E-B), 



(9b) 
(9c) 



e e (x) 



\eB\ 
2nh 

E + E F 
2ixhw c 



m 



- X) dx 
X 



arcsm ■ 



x-X 



— sin I 2 arcsin ■ 
2 V R c 



(9d) 
(9e) 



9 h = Q e (-E,-B). 
Here the cyclotron radii i? 0j h and the classical turning 
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FIG. 2: Classification of the possible orbits. The blue/red lines correspond to electron/hole trajectory segments. The black 
arrows show the direction of the group velocity (cf. the slope of the dispersion curves in Figs. [3] and [4| for a magnetic field 
pointing out of the plane. 



points T± h for electrons and holes are given by 

i? c ,h = Roy/l ± E/E F , R c = k F l 2 , (9f) 

t% = mm{W,X + R e }, (9g) 

T e _ = max{0,A - R e } , (9h) 

r h ± = tU-B,-E), (9i) 

where X — sign(e_B)/cZ 2 is the guiding center coordi- 
nate. Finally, contributions to the Maslov index of a 
given type of orbit come from the collisions with the wall 
of the waveguide, from the collisions with the supercon- 
ductor (Andreev reflections), and from the caustics of the 
electron (hole) segments of the orbit (see Tablc|l|. 

IV. DISPERSION RELATION AND THE PHASE 
DIAGRAM FOR NS SYSTEMS 




k/k F 



FIG. 3: The energy spectrum obtained from Eqs. (|8]) 
(solid/blue line) and (0) (red crosses) as a function of k. The 
parameters are fo = 40, Ao/ft^c = 2, kpW = 106.7 and 
Rc/W = 0.375. In the semiclassical calculations only the 
orbits of type A, E and G need to be taken into account. 



In this section we compare the numerical results ob- 
tained from the exact and from the semiclassical calcula- 
tions outlined in Sections HT1 and ITTT1 The above discussed 
classical orbits can exist in different parameter ranges, 
depending on the strength of the magnetic field, on the 
Fermi energy and on the width of the normal lead. It 
is convenient to use the following dimensionless param- 
eters: i/q, A /Huj c , k^W and R c /W . Figures [3] and 0] 
show comparisons of the exact quantum calculations with 
the semiclassical results for the dispersion relation of the 
quasiparticles. In case of Fig. [3J the magnetic field is 
strong so that the cyclotron radius is smaller than the 
width of the waveguide. One expects therefore that type 
A orbits and Landau levels (corresponding to orbits of 
type E and G) would appear in the spectrum. One can 
see that this is exactly the case, the Landau levels ap- 
pearing as dispersionless states. The agreement between 
the quantum and semiclassical calculations is excellent 
except in the transition regime between type A and type 
E (G) orbits. In case of Fig.[J]thc magnetic field is weaker 
than for Fig. [3] and therefore R c is now larger than W . 
No Landau levels appear and the dispersion relation can 
be well approximated semiclassically using orbits of type 
B, C, D, F and H. 

Whether or not a given type of classical orbit is allowed 
for a certain set of parameter values depends on the po- 
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FIG. 4: The same as in Fig. [3] with parameters i/o — 160, 
Ao/hu c = 8, k F W = 106.7 and R c /W = 1.5. In the semi- 
classical calculations the orbits of type B, C, D, F and H are 
involved. 



sitions of the turning points with regard to the wall of 
the wave guide and to the NS interface. The conditions 
for each type of orbits are summarised in Table Q] The 
turning points depend [see Eqs. j9f]) - ((9l|] on the width 
of the lead W, on the Fermi wavenumber fcp, on the mag- 
netic field (or equivalently, on the cyclotron radius R c ), 
and on the wave number k. In an experiment the for- 
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type of orbit 


N(E) 




conditions for orbits 


A 


S e (rl,rt)-S h (4 


rl) 


iarccos(f ) 






B 


S e (rt,Tt) - S h (4 


rl) 


-§ + Iarccos(§) 


A 


= W, t ] 1<W, rt,r!t=0 


C 


S e (rl,Tt) - S h (4 


rl) 


-i + ^arccos(f) 




<W, tX = W, t1,t!1=0 


D 


S e (rl,rl) - S h (4 




Iarccos(f ) 




tI,t\ = W, t±,t!L = 


E 






l 

2 




T%<W, Tt > 


F 


S*(t',t1) 




1 

4 




= W, rl>0 


G 


~S h (4,r1) 




1 

2 




t\ < W, rh. > 


H 


-S*(t£t-) 




3 
4 




t+ = W, > 



TABLE I: N(E) related to the actions defined in Eq. (9), and the conditions for the possible orbits shown in Fig [2] 



mer two parameters would be fixed, R c could be varied 
by varying the magnetic field, while the wave number k 
of the injected electrons would be uniformly distributed 
if there are many open channels in the quantum point 
contact. For a given W and fcp, ie for a given experi- 
mental sample, one can then draw a "phase diagram", 
which shows the allowed types of classical orbits as a 
function of R c and k. An example of such a diagram is 
shown in Fig. [5l Note that the weak energy dependence 
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FIG. 5: The phase diagram of the allowed type of orbits 
as a function of R c and k. The white regions are classically 
forbidden. The energy of the quasiparticles is E = 0. 

of the turning points translates into a similarly weak en- 
ergy dependence of the phase diagram. Thus for Fig. 
we have chosen E = 0. The white regions are classically 
forbidden, as no classical orbit can be realized for these 
parameter values. 



V. COHERENT ELECTRON-HOLE FOCUSING 

Having obtained the spectrum of the quasiparticles, we 
have now all the necessary information to calculate the 
conductance between the injector and collector for the 
NS system shown in Fig.Q] For Wi, W c <C X F the calcu- 
lations in principle could be carried out using the exact 
results of Section [TTJ However, this would involve eval- 
uations of the Whittaker functions [see Eq. which 
would render the numerical calculations rather difficult. 
Therefore we calculate the conductance semiclassically, 



adopting and generalising the method of Ref. Il2| to ac- 
count for all types of current carrying modes. Namely, 
the dynamics of hole type quasiparticles created by An- 
dreev reflections also needs to be taken into considera- 
tion. We assume that Wi , W c <C R c meaning that the 
angular distribution of the injected electrons is not per- 
turbed by the magnetic field. 

In classical picture injected electrons having orbits of 
type B or D can contribute to the conductance, since only 
these orbits can reach the collector (note that for type F 
orbit the group velocity points in the — y direction). As- 
suming that the wave function in the wave guide is un- 
perturbed by the presence of the collector—, the current 
at the collector is given by 



<9* c 



dx 



(W,L) 



5* h 



dx 



(W,L) 



(10) 



where e is an yet undetermined parameter but will drop 
out when we calculate the conductance. This expression 
is the generalisation of Eq. (16) in Ref. by including 
the contribution of the holes. In WKB approximation the 
wave function ^^(x,y) = (^ c (x,y),^fh(x,y)) T in the 
wave guide is the sum over all classical trajectories from 
the injector to the point (a;, y) of an amplitude factor 
times a phase factor. As in Ref. [Tl we transform the 
sum over trajectories into a sum over modes using saddle 
point integration. Finally, we find 




A* n e lk * L -™ cos af, (11a) 



K=B,D 



and 



9*h 

dx 



(W,L) 



A D & ik n L-i-K cos ^.(1 lb) 



For simplicity, here we give the definitions of the different 
terms appearing in the above expressions only for elec- 
trons having type B orbits. For holes and for type D 



orbits similar expressions were derived but are not pre- 
sented here. The wave numbers k n of the excited modes 
(for a given magnetic field) can be obtained from the dis- 
persion relation by solving the equation E(k n ,R c ) = E. 
The amplitudes A^ n of the modes related to type B of 
orbits are given by 



At 



d(a r , 



' L cos a n 



2v F L d'(a)\ a = 



(12) 



where Ii is the current injected from the injector, v~f = 
kkp/m is the Fermi velocity of the quasiparticles and 
the prime denotes derivation with respect to a. Here a, 
a n and d{a n ) are defined in the following way: the dis- 
tance between two subsequent bouncing at the wall of the 
wave guide for an electron injected at angle a (measured 



from the y axis) is d(a) = 4yi? 2 — (W — i? c sina) — 
2R C cos a. Then p = L/d{a) is the number of bounces be- 
tween the injector and collector. The angle a can also be 
expressed by sin a = (W — X)/R c and it can be related to 
the wave number k n of the modes. Since for mode k n the 
guiding centre coordinate is X = sign(e-B)fc„/ 2 , we have 
sina„ = W/R c — sign(eB)k n /kF (because R c = k-pl 2 ), 
and then p n = L/d(a n ). Finally, Sf ^(p n ) is related to 
the action of electrons (holes) calculated from the injector 
to the collector. In case of type B orbits the summation 
over n includes only those modes for which the group ve- 
locity dE(k) / d(hk) is positive, i.e the mode propagates 
from the injector to the collector. (Note that the group 
velocity is basically given by the slope of the curves in 
FiggJ. 

The conductance G(E, R c ) between the injector and 
the collector is G = Ii/V c , where V c is the collector volt- 
age. Taking into account Eq. (fTO]) it can be calculated 12 

as 



V, = e 



Go 



\ 9 -t(W,L)\ 



(13) 



Here G e (Gh) is the conductance of the point contact 
of the collector for electrons (holes). It is estimated by 
Eq. (17) in Ref. [ID- Since both G e and Gh are propor- 
tional to the parameter e, it drops out from V c . Similarly, 
the injected current 7j drops out from G(E, R c ) since the 
derivatives of the wave functions in (|lip are proportional 
to \JT~i through the amplitudes A Pn . 

Equations (jTUJ) , (fTTj) and (fT3"| allow us to calculate 
the conductance semiclassically as a function of the en- 
ergy of quasiparticles and of the magnetic field. This 
is the main result of our paper. An example is shown 
in Fig. [SI One can clearly see that for small magnetic 
fields (W/R c < 0.4) the conductance is a rapidly oscil- 
lating function of the magnetic field and can also be neg- 
ative. The first observation is the consequence of the con- 
structive interference of many coherently excited modes 
(edge states) . Similar effect was found by van Houten et 
all 12 ' 13 . The fact that the conductance can be negative 




FIG. 6: The conductance G for low magnetic fields at E = 0. 
Inset: the conductance for stronger magnetic fields. We used 
k F W = 26.7. 



is a consequence of the hole-like excitations in the hybrid 
NS system. Indeed, at low magnetic fields (large R c ) an 
injected electron will undergo one or more Andreev re- 
flections and it can happen (see type D orbits) that a hole 
will arrive to the collector resulting in negative conduc- 
tance. This is the so called Andreev-drag effect^ . The 
heights of the positive peaks are comparable with those 
of the negative ones, indicating that it can be a robust 
effect. This observation is in line with the findings of 
Ref.0 and also with the results of Ref. 17 where an exact 
(numerical) quantum calculation has been performed for 
the same system. 

The use of WKB approximation means that our re- 
sults should be accurate at low magnetic fields when a 
large number of edge states are populated. Nevertheless, 
we find that this method, at least qualitatively, also de- 
scribes the regime of high magnetic fields. Namely, from 
the inset of Fig. [6] one can see that increasing the mag- 
netic field the conductance decreases and rapidly goes 
to zero for W/R c > 1.5. This can be understood from 
classical considerations. Increasing the magnetic field the 
cyclotron radius decreases and at a certain value of the 
field we have 2R C < W meaning that no injected elec- 
tron can hit the superconductor and undergo Andreev 
reflections. Instead, the electrons move to the left skip- 
ping along the wall of the wave guide (type F of orbits) 
and eventually they leave the system without reaching 
the collector, ie the conductance becomes zero. Accord- 
ing to our semiclassical calculations, for a system with 
k-pW = 26.7 as in case of FigJS]the last current-carrying 
mode disappears when W/ R c > 1.8, in broad agreement 
with the classical picture. 



VI. CONCLUSIONS 

In conclusion, we have studied the transverse electron- 
hole focusing effect in a normal-superconductor system 
similar to the setup of Ref. d. As a first step to determine 
the conductance, we calculated the energies of the bound 
states both quantum mechanically and in semiclassical 
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approximation. We have shown that semiclassical meth- 
ods well reproduce the results of the relevant quantum 
calculations and thus can help to understand the under- 
lying physics. We have identified those classical orbits 
which contribute to the conductance in our system and 
derived a semiclassical formula for the conductance. In 
agreement with the quantum calculations of Ref. for 
weak magnetic fields the semiclassical conductance shows 
rapid oscillations and the presence of the Andreev-drag 
effect. For stronger magnetic fields the conductance goes 
to zero which can be understand invoking the classical 
dynamics of electrons at such fields. Our results can be 



considered generalisation of similar work o 12 ' 13 for normal 
systems since in our case the current carrying modes are 
comprised of electron-like as well as hole-like quasiparti- 
cles. 
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